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Cluster number counts can constrain the properties of dark energy if and only if the evolution in the 
relationship between observable quantities and the cluster mass can be calibrated. Next generation 
surveys with ~ 10 4 clusters will have sufficient statistics to enable some degree of self-calibration. 
The excess variance of counts due to the clustering of clusters provides such an opportunity and can 
be measured from the survey without additional observational cost. It can minimize the degradation 
in dark energy constraints due to an unknown power law evolution in the mass-observable relation 
improving constraints on the dark energy equation of state by a factor of 2 or more to o~(w) = 0.06 
for a deep 4000 deg 2 survey. 



I. INTRODUCTION 

The abundance of clusters of galaxies of a known mass 
is a sensitive probe of linear density fluctuations given 
Gaussian random initial conditions (e.g. Q). Since the 
growth of fluctuations is halted by the dark energy, abun- 
dance measurements as a function of redshift can be 
translated into dark energy constraints 0, Q ■ 

Unfortunately, the mass of a cluster is not a direct ob- 
servable. The key to exploiting this dark energy sensitiv- 
ity is to determine the relationships between observable 
quantities such as flux or temperature and the mass. Ab 
initio computations from cosmological simulations serve 
as useful guides but invite misinterpretation when used 
directly due to missing gas, star formation, and AGN 
physics (e.g. 0,0)- Much like with the distance ladder 
determinations of the Hubble constant, more progress 
can be made by cross-calibrating the mass-observable re- 
lations between X-ray, lensing, optical and microwave 
surveys. A potential drawback is that cross-calibration 
makes the interpretation of each data set subject to the 
systematic errors in all of the surveys employed and so 
the weakest link in the chain. Cross-calibration can also 
only be performed across the redshift and mass range 
where surveys overlap. 

Recently, self-calibration techniques have been advo- 
cated as a useful check on cross-calibration and the in- 
ternal consistency of dark energy determinations. Clus- 
ters benefit from having several mass-sensitive proper- 
ties that can be reliably predicted from simulations of 
their dark matter properties only. The first of course 
is the evolution of the abundance itself. An unknown 
constant normalization of the mass-observable relation 
can be determined from the evolution in the abundance 
above threshold in the given observable jfj. Even an 
arbitrarily evolving normalization can be determined by 
measuring the abundance as a function of the threshold 
and redshift 0. This technique however requires a dy- 



namic range in the mass function that is substantially 
larger than the scatter in the mass-observable relation 
and will only be possible with very deep surveys or at 
fairly low redshifts. Finally, mild evolution in the mass- 
observable relation can be calibrated by examining the 
power spectrum of the clusters |8( ■ The drawback of the 
latter method is that measuring the three dimensional 
power spectrum generally requires precise redshifts and 
hence costly spectroscopy or an improvement in photo- 
metric redshift techniques beyond that required for abun- 
dance studies themselves (Az ~ 0.1). 

Here we point out that there is an analog of the power 
spectrum self-calibration technique that will automati- 
cally come out of the statistical analysis of the number 
counts in any given cluster survey. The excess or sample 
variance of the counts due to the clustering of clusters 
must be included in future error analyses [9j- For self- 
calibration purposes, this source of "noise" is actually 
signal and can be used to calibrate the mass observable 
relation. 



II. NOISE AS SIGNAL 

The probability of measuring a number of clusters N 
in a cell of a given redshift z and angular extent is given 
by the Poisson distribution 



P(N\m) 
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where the mean m = (N) p and the brackets denote av- 
eraging over realizations of the Poisson process. 

Now take a set of cells indexed by i where the mean 
number m,i fluctuates in space 



d 3 xWi(x)n(x; Zi 
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where Wi(x) is the cell window function, in this case a 
top hat, and n(x; Zi) is the spatial number density. On 



2 



large scales, fluctuations in their spatial number density 
trace the linear density fluctuations <5(x; z) from the large 
scale structure of the universe 

n{x;z)=n(z)[l + b(z)5(x;z)], (3) 

where b is the linear bias of the clusters. Overbars denote 
a spatial average or a sample average over realizations 
of the large scale structure. Thus the sample averaged 
number counts rfii = Vifi(zi) where the cell volume is 
Vi = J d?xWi. 

The mean numbers then possess a sample covariance 
given by the linear power spectrum P(k) 

S i;j = ((m,i - fhi)(rrij - rhj))s 
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Here W^k) is the Fourier transform of Wj(x) and 6j = 
b(zi). For a single spherical cell of comoving radius R, 

1 /2 

the fractional errors /rrii = bi<jR where an is the rms 
linear density fluctuation in the cell. 

The likelihood of drawing a set of cluster counts N = 
(Ni, ...N c ) given a model for rft and S is then 



L(N|m,S) = J 



d c m 



i=l 



G(m|m,S), (5) 



where G denotes the multivariate Gaussian distribution 



G(m|m,S) 



1 



tj 1 (m— m)S 1 (m— m) 



(6) 



v/(27r) c detS 

It is instructive to consider a few special cases. In the 
limit that sample variance is negligible compared with 
the Poisson variance Sa <C rhi the Gaussian approxi- 
mates a delta function and 



L(N|m,S)«nP(7V J |m l ). 
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This form is appropriate for very rare clusters and is used 
in the analysis of local high temperature clusters. 

In the limit of large numbers rrii 3> 1, the Poisson 
distribution approaches a Gaussian 



][[F(iV l |m l )«G(N|m,M), 
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where M = diag(nij). The likelihood becomes a convo- 
lution of Gaussians 



L(N|m,S) 



d c m G(N|m, M)G(m|m, S) 
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or via the convolution theorem approximately a Gaussian 
inN-m with covariance C = S + M. 

The statistical properties of the counts are hence speci- 
fied by their mean m and their sample covariance S. The 
space density n(z) controls m and hence the "signal" ; the 
bias b(z) controls S and hence the "noise" . Given a cos- 
mology, both n and b can be predicted as a function of 
the cluster mass (see fllVf) . From the perspective of mass 
calibration both ingredients are therefore signal. 
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FIG. 1: Mass sensitivity of the counts and their variance compared 
with sensitivity to the dark energy equation of state w at z = 0. 
The integral of the mass function dn/dlnM, here normalized to 
Mo = 10 14 /i — 1 Mq from the threshold Aftb predicts the expected 
counts. Changes due to w can be compensated by a comparable 
shift in M t h (dark vs. light shaded region) but also change the 
fractional cell variance at threshold (horizontal lines), shown here 
normalized to the typical cell volume (50h _1 Mpc radius) as fe 2 cr| . 



III. FISHER MATRIX 

The Fisher matrix formalism provides a means of mak- 
ing projections for the result of likelihood analyses of fu- 
ture data given a parameterized model for the mean m 
and sample covariance S of the counts as well as a fiducial 
choice of the true parameters. 

Given a set of parameters p a on which these quanti- 
ties depend, the information is quantified by the Fisher 
matrix 



_ dHnL 
f a p = — ^ — / 

OPaOpp 
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such that the parameter covariance matrix C a p w 
(F _1 ) Q/ 3. The marginalized parameter errors are then 

a(p a ) = Caa. 

Again let us begin by considering the limiting cases. 
In the case of negligible sample covariance 

F a p = —fhi^rhifi = m^M^m ^ , (11) 

where , = d/dp a ^(j- ^ n the case of negligible Poisson 
errors 11] 

F a(3 = m^S-'m^ + ^TrfS-^S-^] . (12) 

Notice that even in the rrii ^ 1 Gaussian limit, the Pois- 
son Fisher matrix does not carry a term involving the 
derivatives of the variance M )CB . However the fractional 
error induced by including such a term scales as m~ 1 and 
is therefore negligible in this limit. 
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Given these two limits, we approximate the full Fisher 
matrix as 

F a fj = m^C" 1 ^ + i , &[C- 1 S,aC- 1 S, / j] , (13) 

where recall C = S + M. The two pieces represent 
the contribution to the information on the parameters 
from the mean of the cell counts and their cell-to-cell 
(co)variance induced by structure in the universe. 



IV. SELF-CALIBRATION 

The two sources of information, the mean counts and 
their cell-to-cell variance, depend differently on the clus- 
ter mass. By comparing the two one can in principle 
solve for the mass and hence remove the calibration un- 
certainty in the interpretation of the number counts. 

Given an initial power spectrum, simulations can re- 
liably predict the number density of dark matter halos 
associated with clusters of a given mass. For illustrative 
purposes, we will employ the fitting function 12] 



dn p. 
= 0.3 P 
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cxp[-|lno-" 1 + 0.64| 3 - 82 ] , (14) 



where a 2 (M;z) = a\(z), the density field variance in 
a region enclosing M = 4irR 3 p, n /3 at the mean matter 
density today p™. Likewise the bias of these objects can 
be described bv pLIT^ 



b(M; z) = l + 



aS 2 /a 2 1 



2p 



S c [l + (aS 2 c /a 2 )P] 
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with a = 0.75, p — 0.3, and S c — 1.69. Though these 
relations must be replaced by numerical results in a real 
analysis, self-calibration will still be possible so long as 
the bias and mean number counts scale with the mass of 
the objects in a predictable way. 

Consider now a selection that is defined by a threshold 
in some observable quantity such flux or temperature. 
Let us define the relationship between the observable 
threshold hh and the mass threshold by two parameters 
A and m M 
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where /o is an arbitrary normalization parameter, M (z) 
characterizes an a priori guess for the relation such that 
deviations can be described by a power law in (1 + z), 
and p is considered known (see Q for a generalization to 
arbitrary p(z)). Then the statistical model of the counts 
is defined through Eqn. © and Q 



n(z) = 



rflnM 



dn 



In M th (z) 



dhxM 
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b(z) = - / d l nM -p-b(M;z). (17) 

n J In M th (z) dlnM 



In reality the mass-observable relation will contain finite 
scatter which will blur the threshold; this scatter must 
also be modeled in a real analysis. We use this simple 
prescription for illustrative purposes only. 

Fig. illustrates the self-calibration idea. If only the 
counts are considered, changes in the cosmology are de- 
generate with those in the threshold. However lowering 
the threshold to compensate a smaller amplitude of den- 
sity fluctuations has two effects on the variance of counts 
that break the degeneracy: it lowers the variance due to 
the decreased underlying linear structure but also makes 
objects at a fixed mass rarer and hence more highly bi- 
ased. 

We approximate the results of a joint likelihood anal- 
ysis of the mass-observable and cosmological parameters 
via the Fisher matrix. We take the cosmological param- 
eters as the normalization of the initial curvature spec- 
trum <5 f (= 5.07x 10~ 5 ) at k = 0.05 Mpc^ 1 (see for its 
relationship to the more traditional as normalization), its 
tilt n(= 1), the baryon density H.bh 2 (— 0.024), the dark 
matter density fl m h 2 (= 0.14), and the 2 dark energy pa- 
rameters of interest: its density finE(= 0.73) and equa- 
tion of state w{= —1) which we assume to be constant. 
Values in the fiducial cosmology are given in parentheses. 
The first 4 parameters have already been determined at 
the few to 10% level through the CMB and we will 
extrapolate these constraints into the future with priors 
of cr(ln^) = a{n) = a(lnn b h 2 ) = a(lnfl m h 2 ) = 0.01. 

For illustrative purposes let us take a fiducial clus- 
ter survey with specifications similar to the planned 
South Pole Telescope (SPT) Survey: an area of 4000 
deg 2 and a sensitivity corresponding to a constant Mq 
= Mth|fid = 1O 14 - 2 /i _1 M . We divide the number counts 
into bins of redshift Az = 0.1 and 400 angular cells of 
10 deg 2 and vary the maximum redshift z max for which 
photometric redshifts will be available. With these large 
cell sizes, the covariance between neighboring cells is neg- 
ligible, considerably simplifying the analysis. 

Constraints in the Ode — w dark energy plane from 
counts alone are severely compromised by mass thresh- 
old uncertainties (see Fig. EJ . This is especially true for 
the lower z max since the lever arm is not sufficient to dis- 
tinguish between the power law deviations in the mass- 
observable relation and the dark energy. Here the errors 
in f^DE degrade by factor of 37 and w by a factor of 
4.6. Much of the dark energy information is restored by 
including variance information. For z max = 1, errors im- 
prove by a factor of 7.5 to o^JIde) = 0.04 and a factor of 
2 to a(w) = 0.09. 

Self-calibration also makes the errors more robust to 
Zmax, for z max = 2, a(fluE) = 0.03 and a(w) = 0.06. 
These results are also relatively insensitive to the cell size. 
For 4 deg 2 cells or Az = 0.05 the errors remain nearly 
the same. Likewise, the relative improvement due to self- 
calibration is relatively insensitive to the mass threshold 
assumed. Finally errors on the actual calibration param- 
eters are a(A) = a(m) — 0.14 for z max = 2. 
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FIG. 2: Projected constraints in the Que — w dark energy plane 
(68% CL) for the fiducial survey and a maximum redshift of (a) 
Zmax = 1 and (b) 2 ma x = 2. Unknown power law evolution in the 
mass-observable relation degrades constraints from the inner most 
to outermost ellipses if only the abundance or counts is used. By 
adding in information from the count variance, the survey partially 
self-calibrates leading to the middle ellipses. 



V. DISCUSSION 

Cluster number count surveys contain information not 
only in the mean counts but also in their cell-to-cell vari- 
ance. The latter depends on the clustering of clusters and 
hence provides an independent constraint on their mass. 
This information automatically comes out of a full likeli- 
hood analysis of the counts and provides an opportunity 
for self-calibration of the survey. 

In the idealized case where the scatter in the mass- 
observable relation is low and known, self-calibration 
via the variance improves dark energy constraints for 
power law evolution in the mass threshold by a factor 
of 2-10 depending mainly on the maximum redshift for 
which accurate photometric rcdshifts can be obtained. In 
principle, self-calibration can constrain a more arbitrary 
evolution, e.g. an independent mass-observable relation 
for each redshift slice. In practice, the resulting dark 
energy constraints are then too weak to be of interest 
(<j(Q£,e) — 0.07, a(w) = 0.57 in the fiducial survey). 
Likewise, for an unknown scatter in the mass-observable 
relation, variance self-calibration alone is unlikely to suf- 
fice. In these cases it can be supplemented with selections 
at various thresholds in the observable Q or with the full 
angular power spectra of the cell counts (e.g. though 
a full treatment is beyond the scope of this work. 

Since self-calibration involves only information which 
exists in the survey itself, it comes at no additional obser- 
vational cost. It therefore complements potentially more 
precise but costly cross-calibration studies. 
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